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Abstract — It is already well-known that interference alignment 
(IA) achieves the sum capacity of the if -user interference channel 
at the high interference regime. On the other hand, it is clear 
that in the very low interference regime, K (as opposed to if/2) 
degrees of freedom (DoF) can be accessed by simple ("myopic") 
single-link multiple-input multiple-output (MIMO) techniques 
such as waterfilling. Recent results have indicated that in certain 
low-to-moderate interference cases, treating interference as noise 
may in fact be more preferable. In this paper, we present a 
distributed iterative algorithm for if -user MIMO interference 
networks which attempts to adjust itself to the interference 
regime at hand, in the above sense. The proposed algorithm 
combines the waterfilling solution with joint mean squared error 
minimization to unveil the available DoF of the interference 
channel and to perform IA, as necessary. Sum-rate computer 
simulations for the proposed algorithm over Ricean fading chan- 
nels show that, in the interference-limited regime, our algorithm 
reconfigures itself in order to achieve IA whereas, in the low-to- 
moderate interference regime, it leads itself towards interference- 
myopic MIMO transmission. 

I. Introduction 

Interference alignment (IA) is a recently proposed transmis- 
sion technique for the if -user interference channel which is 
shown to achieve a sum-rate multiplexing gain of if/2 at the 
high interference regime (TJ, J2). IA is based on appropriate 
linear precoding at the transmitters, aiming at post-receiver 
processing interference cancellation, and requires only global 
channel state information at all participating transceivers. 

Exploiting the space dimension of multiple-input multiple- 
output (MIMO) systems to perform IA, several research 
works present IA-achieving precoding designs [1|-|10| and 
investigate the feasibility of IA IfTTI . Ifl2l for the if -user 
MIMO interference channel. For the special case of K = 3, 
a closed-form solution for IA was presented in [2] that was 
further processed in [10] for increased sum-rate performance. 
However, for K > 3 MIMO communicating pairs, closed- 
form solutions for IA are in general unknown and several 
iterative algorithms have been recently proposed (see e.g. 
[3 1 — [ 1 1 and references therein). The vast majority of those 
algorithms targets at implicitly achieving IA through the 
optimization of a constrained objective function. To this end, 
several objective functions have been considered, such as for 
example: f) minimization of the total interference leakage @, 
H], if) minimization of the sum of least squared errors 1131 . 
Hi) minimization of the mean squared error (MSE) Q, (6), 
[|9l , iv) maximization of the signal-to-interference-plus-noise 



ratio (SINR) J5), and v) maximization of the sum-rate 
performance ffl, |7l, iflOl 

Although IA attains the optimum rate scaling at the high 
interference regime, there are certain combinations of SINR 
regimes and channel conditions where it does not [14|-[16|. 
For example, the authors in lfT31 analyzed conditions for 
the intended and interference MIMO channels under which 
treating interference as noise at the receivers is sum-capacity 
achieving. Very recently, the sum-rate performance results of 
lTT6l . for if -user MIMO networks with asymmetric average 
powers and line-of-sight (LOS) conditions among the intended 
and interference channels, demonstrated certain regimes where 
interference-myopic MIMO transmission yields superior per- 
formance to IA. 

Based on the above, it would be desirable to devise a 
sum-capacity-achieving transmission design for the if -user 
MIMO interference channel that is transparent to the interfer- 
ence conditions. To the best of our knowledge, the available 
transmission techniques for the if -user interference channel 
must know a priori the interference conditions so as to choose 
between the two extremes: treating interference as noise or 
performing IA. Inspired by the distributed IA algorithm of 1 5 1 
and the results of 1161 . in this paper we present a distributed 
iterative interference handling algorithm that combines the wa- 
terfilling (WF) solution [17 1 with the minimum MSE (MMSE) 
criterion to unveil the available degrees of freedom (DoF) 
of the interference channel and to perform IA, as necessary. 
Early numerical evidence corroborates our expectation that 
the proposed technique reconfigures itself so as to allow the 
maximum attainable sum rate for the interference channel 
conditions at hand. 

Notation: Vectors and matrices are denoted by boldface 
lowercase and capital letters, respectively. The transpose con- 
jugate and the determinant of matrix A are denoted by 
A H and dct (A), respectively, whereas [A],.j represents the 
element of A. I n is the n x n identity matrix and 
diag{a} represents a diagonal matrix with vector a in its main 
diagonal. The expectation operator is denoted as E{-} whereas, 
X ~ CAf (/i, er 2 ) represents a random variable following the 
complex normal distribution with mean /i and variance a 2 . 

II. System and Channel Model 

We present below the if -user MIMO interference system 
model and the wireless channel model under consideration. 



A. System Model 

A multiuser MIMO system consisting of A pairs of com- 
municating users is considered. In particular, each transmitting 
user (Tx) k, where fc = 1,2, . . . , A, equipped with n!?' an- 
tennas wishes to communicate with the -antenna receiving 
user (Rx) k. All A simultaneous transmissions of symbols 
Sfe G C dfcXl Vfc, with dk < min(n^ , rip ), are assumed 
perfectly synchronized and each Tx k processes with a 
linear precoding matrix G C"t xdk before transmission. 
For the transmitted power per Tx k E{||VfcP^Sfe|| 2 } < P 
with P being the total power constraint per Tx and = 
diag{[P 1 (fe) P 2 (fc) . . . Pjfc*]}, where n = 1, 2, . . . , d k , denotes 
the power allocated to the nth data stream at Tx k. Without 
loss of generality, it is assumed throughout this work that 
E{sfeS^} = Id k Vfc. The baseband received signal at Rx k 
can be mathematically expressed as 

K 

y k =H fe!fe V fc P|s fe + H M V^P|sf + n fe (1) 



where H M G C n R xri T, with i = 1,2,..., K, denotes 

the channel matrix between the fcth Rx and the l\h Tx, and 

[fe] 1 

rifc G C™r represents the zero-mean complex additive white 
Gaussian noise vector with covariance matrix ail [ k \. After 
signal reception, each Rx k is assumed to process y^ with a 
linear filter U fe G C n R xrffc as U¥y k . 



B. Channel Model 

The flat fading channel model of lfl6l is assumed for which 
the channel matrix between Tx I and Rx fc is given by 



H 



H Mj _ fc = t 



(2) 



where parameter a kl i G [0, 1] is used for modeling asymmetric 
average powers among the intended and the interference chan- 



nels, and Hk,e G 
defined as 



1 t , which describes Ricean fading, is 



H 



i 



KkJ + 1 



HNLOS 



In (O, K k ,e is the Ricean K-factor and H^ os G 
is the scattered component of Hk,e sucn that [H^^ s ] t . 3 



CAf(0,l) Vi = l,2,...,n^ ] and Vj = l,2,...,n^ J . More- Vfc = 1,2. 
over, a e (6 t ) G C n T xl and a k {9 r ) G C"r x1 denote the 
specular array responses at Tx I and Rx fc, respectively, with Q t 
and 9 r being the angles of departure and arrival, respectively. 

III. An Interference-Adjustable Distributed 
Algorithm for A'-user MIMO Networks 

In this section we provide the motivation and formulation of 
the proposed distributed iterative algorithm. A brief discussion 
of the algorithm's convergence is also included in the end. 



A. Motivation 

As mentioned earlier, the sum-capacity-achieving transmis- 
sion techniques for the Abuser MIMO interference chan- 
nel strongly depend on the interference conditions. Treating 
interference as noise is better in weak interference levels 
[14|-[16|, whereas IA achieves the sum capacity at the high 
interference regime (TJ, 0. To this end, choosing between 
the latter two transmission techniques requires the a priori 
knowledge of the interference conditions. Furthermore, all 
IA-achieving algorithms require IA feasibility conditions to 
be met a priori (see e.g. JT] — [ 10|). For example, for K 
simultaneous riR x nx MIMO transmissions, each Tx fc is 
restricted to transmit dk < (riR + n^)/(K + 1) data streams. 
However, [ 16 1 demonstrated low-to-moderate interference sce- 
narios where interference-myopic MIMO transmissions, each 
aiming at the individual user rate maximization, yield higher 
sum-rate performance than IA. For those scenarios, there was 
at least one Tx transmitting more than (nR+nx) / (if +1) data 
streams. Finally, a typical feature of IA-achieving techniques 
is the equal-power allocation on each Tx's data streams |Q]- 
J8). Inspired by the capacity-optimal strategy for single-user 
MIMO systems [ 17] and the sum-rate results for Ricean fading 
channels in ifTBI . we intuitively expect that the equal-power 
allocation per Tx will be suboptimal in the weak interference 
regime and under strong LOS conditions. 

Motivated by all the above, we formulate in the following 
a distributed iterative algorithm that implicitly chooses dk for 
each Tx fc accordingly to the SINR and LOS conditions and 
jointly designs Vfe's, Pfe's and Ufc's for all transceivers in the 
network for increased sum-rate performance. 

B. Algorithmic Formulation 

After processing y k at each Rx fc with Ufe, the block 
MSE for the considered A-user MIMO interference network 
is expressed as 



K 

JMSE = ^E{||UHy fc -s fe || 2 }. 



(4) 



k=l 



In this work we are interested in the joint design of V^'s 
and Ufc's that minimizes (0). Therefore, the considered joint 
MMSE constrained optimization problem is mathematically 
expressed as 



[f] 



min Jmse 
v fc ,u fc 

s.t. ||V fe „|| 2 = ||U fen | 
, K and Vn = 1,2... 



(5) 

2 = 1 

dk- In ©, V fe „ and 
Ujfc„ denote the nth column of Vfc and Ufc, respectively. 
Similar to 0, 0, 0, and assuming reciprocal forward 
and reverse networks, we adopt alternating minimization to 
iteratively obtain Vj,'s and U^'s satisfying (0, as follows. 

1) Forward Network: In the original forward network, 
given Vfc and Pfc V k = 1, 2 . . . , K, the nth column of Ufc at 
Rx fc satisfying (|5]l is obtained as 



B fe 1 Hfe j feVfe r , 
|B., 1 Hfe : feVfc r , 



(6) 



where £ 



M [fc] 

n K xn R j s gj ven by 
A" 



B k = £ H M V,P,V £ H H^ + ofr^m . 



2) Reciprocal Network: In the reciprocal network, each Rx 
k utilizes Ufe to transmit to its intended Tx k 0. Then, each 
Tx k computes and P k through the following two-step 
procedure. 

Step I: Given U fe and P k Vfc = 1,2..., if, each Tx k 
computes its receive filter Efe G C™ T x k so as to satisfy a 
reformulated version of © for the reciprocal network. Hence, 
the nth column, with n = 1, 2, . . . , dk, of Efc satisfying © is 
derived as 



E 



(B/c) Hfc.fcUfer, 



Algorithm 1 Interference-Adjustable Joint MMSE 

1: initialization: Set dk = min(n!^ , ) Vfc = 1,2,...,^ 
(7) and start with arbitary precoding matrices £ C" T xdfc 

and power allocation matrices P k = ~P/dkJ- n [k] 
2: Begin iteration 

Forward Network 
3: Compute Bfc at each Rx k according to © 
4: Obtain each Rx k MMSE-based filter U k using © 
Reciprocal Network 
Step I: 

Compute Bfc at each Tx k according to © 
Obtain all Tx MMSE-based filters E fc 's using © 
Step II: 

Compute Cfc at each Rx fc according to (TTOb 



kn 



where B&G C™t x ™t is given by 



Bfe) HP,Ufe„ 



(8) 



Bfe= HfjJJiPtUf U itk + crfel n w 



(9) 



5: 
6: 

7: 
8 
9 
10 

II: 



Perform SVD to each C^Hfc^Efc according to (ITTb 
Obtain each Tx k precoding matrix ~V k using ( fT2l 
Compute Pj. for each Tx fc from the WF solution for 
channel C k 5 H fcjfc E fc 

Repeat until convergence, or until the number of iterations 
reaches a predefined limit 



Step II: After obtaining all ~E k 's from © and given P? 
V I ^ k, each Tx k computes its optimum precoding and power 
allocation matrices for its effective channel Hfc^.Efc, assuming 
knowledge of the interference plus noise covariance matrix 
Cfc G C™r "r , which is given by 



A" 

E 

= l,t^k 



Hfe j»E£PfE?H H 



fc,£ 



k r 



(10) 



In particular, the optimum precoding matrix F k G <C dkXdk for 

the MIMO channel H k . k ~E k knowing C/j is obtained from the 

_ i 

singular value decomposition (SVD) of C fe 2 H k . k ~E k , i-e- as 



WfeAfcF" 



(11) 



"rX^r 1 and A fc G C"r xdfc . Then, using F fc 



where g 

and the MMSE-based filter Efe, Vfe for Tx k is obtained as 



V fc =E fc F fc . 



(12) 



Finally, the power allocation matrix Pfe at each Tx k is derived 
from the WF solution ifTTl for the channel C fe 2 H k _ k ~E k . 

The proposed distributed iterative algorithm for the if -user 
MIMO interference channel is summarized in Algorithm Q] 

C. Convergence 

In each iteration the proposed algorithm attempts to min- 
imize the global objective function ,7mse via alternating 
minimization. Similar to several alternating minimization tech- 
niques ©, ©, (9|, our numerous computer simulation results 
showed a per iteration MSE reduction and we often observed 
global convergence. A thorough convergence analysis of the 
proposed algorithm will be the objective of future work. 

The proposed algorithm reconfigures itself to the interfer- 
ence and channel conditions to effectively choose d k and 



design and P^ at each Tx k. This is accomplished via 
the WF solution adopted in Step II. In particular, the number 
of data streams per Tx k is initialized as dk = min(n!p',np'). 
Then, at each iteration, the subspace minimizing Jmse at each 
Tx is computed. Within each such subspace, and P k that 
maximize the instantaneous rate for Tx k are obtained. It is 
noted that, thanks to the following property ifTUl 



span(H feife E fe ) = span(H fcifc V fe P fc ), 



(13) 



the optimum and P k for C fe ^t^Et minimize Jmse- 
Note also that allocating zero power to some streams, i.e. 
reconfiguring the number of data streams d k , seems not to 
influence the per iteration global MSE reduction. 

In order to check the convergence of our algorithm, in each 
iteration we compute the instantaneous achievable sum rate, 
which is defined as 

A rf, / p W n H hjjHu \ 

/,■„• 2l E 1 ' 



fe=l n=l 



where G 



is given by 



A 



E E p i 



HttH 



k.e 



= l,t=£k i=l 



(14) 



(15) 



E 



Note that obtaining the sum rate using (fl4l implies that no 
further receive processing, such as e.g. MMSE successive 
interference cancellation, apart from XJ k is utilized at each Rx 
k. The algorithmic iterations terminate whenever R a converges 
to a specific value with a prespecified accuracy. 
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Fig. 1. Sum-rate performance, Ft, versus transmit power per Tx, P, for 
3-user 4x4 MIMO networks over Rayleigh and Ricean fading with a = 1. 



Fig. 2. Sum-rate performance, Ft, versus transmit power per Tx, P, for 
3-user 4x4 MIMO networks over Rayleigh fading with a = 10 -3 . 



IV. Sum-Rate Performance Results 

This section presents numerical simulations for the sum- 
rate performance of the proposed interference-adjustable dis- 
tributed algorithm for a 3-user 4x4 MIMO interference 
network. For comparison purposes, sum-rate computer simu- 
lations for the IA-achieving algorithms: i) optimized subspace 
IA precoding [10, Sec. IV.B.l] and ii) per stream maximum 
SINR 12 are also shown. In particular, we have simulated the 
ergodic sum-rate performance defined as [6 | 



I? 




det ( I_ w + H M VfcP fc VM fc Rfc 



(16) 

where £ {•} denotes the expectation over all channel realiza- 
tions H M Vfc,£ = 1,2,..., K and R fe e C n R xn R is the 
interference plus noise covariance matrix which is straight- 
forwardly obtained from ( [Tol l by just replacing E/s with the 
V/s of (1121 1. The averaging in ([Tol l was evaluated via Monte 
Carlo simulations for 100 independent channel realizations and 
the channels were normalized as E|||Hfcj-||p| — rJ^nlfi 
Vfc,j = 1,2,..., A". Without loss of generality, for the 
channel model in Sec IH-BI we have assumed that a^j = a 
and Kk,t = K Vfc,l. For all algorithms, Vfc's and UVs were 
randomly initialized with unit norm columns. The proposed 
algorithm was also initialized with = P/4Lt Vfc = 1,2 
and 3 whereas, for both IA-achieving algorithms, IA feasibility 
conditions were set a priori, i.e. dk = 2 V fc = 1, 2 and 3. A 
maximum of 2000 iterations was used per iterative algorithm 
and each of those algorithms was declared converged when 
the difference in its objective function between two successive 
iterations was less than 10 -7 . 

As shown in Fig. [T] for various Ricean fading channels 
with a = 1, the rate of the proposed algorithm is similar 
to that of the per stream SINR maximization. In particular, 



in the noise-limited regime, our algorithm computes Tx filters 
that maximize signal powers over the noise power. The latter 
maximizations include optimum power allocation, resulting in 
a slightly higher R than that of the per stream SINR maximiza- 
tion under weak interference conditions. For Rayleigh fading, 
this behavior happens for P values ranging from — 5 to 5 dB 
whereas, for Ricean fading with k = 10, the gains from power 
allocation are higher and for a wider range of P. To this end, 
optimum power allocation seems to result in higher R gains 
as LOS conditions among all users in the network become 
stronger. More importantly, in the interference-limited regime, 
the proposed algorithm reconfigures itself so as to achieve 
IA, which is rate-scaling optimal under strong interference 
conditions. Figs. [2] and [3] depict R versus P for Ricean fading 
channels with k = and k = 10, respectively, and a = 10~ 3 . 
As expected and clearly depicted in both figures, the IA algo- 
rithms are restricted to achieve a multiplexing gain of K/2. 
On the contrary, the proposed algorithm achieves all available 
K DoF and outperforms both IA-achieving algorithms in the 
entire range of P. It is noted that for this regime the behavior 
of the proposed algorithm resembles that of the interference- 
myopic transmissions with dk = 4 Vfc presented in |16|. 

In Fig. 2] the convergence of R a for Rayleigh fading 
channels with a = 1 and for various values of P is depicted. 
As shown, R a converges fast to a specific value and the speed 
of convergence depends on P. In particular, as P increases, 
more algorithmic iterations are needed for the convergence 
of R a . Similar trends for the convergence of R a have been 
found for Ricean fading channels and various values of a. 
Specifically, additional convergence results have shown that P 
and a influence the i? a 's speed of convergence. As P — > and 
a — > less algorithmic iterations are needed for convergence. 
Moreover, the convergence of R a seems to be independent 
of the Ricean K-factor. It is finally noted that the sum- 
rate simulations in Figs. QJ[3] showed that the convergence 
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Fig. 3. Sum-rate performance, R, versus transmit power per Tx, P, for 3-user 
4x4 MIMO networks over Ricean fading with k = 10 and a = 10 -3 . 



Fig. 4. Convergence of the instantaneous achievable sum rate, R a , for 3-user 
4x4 MIMO networks over Rayleigh fading with a = 1. 



of the proposed algorithm is much faster than that of the 
per stream SINR maximization. For example, in Fig. Q] for 
Rayleigh fading and for P = 30 dB, the proposed algorithm 
achieves IA after 400 iterations whereas, the per stream SINR 
maximization has not converged to the IA solution even after 
2000 iterations. 

V. Conclusions 

In this paper, a novel interference-adjustable distributed iter- 
ative algorithm for the J\"-user MIMO interference channel is 
presented. The proposed algorithm combines the WF solution, 
to unveil the available DoF of the interference channel, and 
iterative joint MMSE, to perform IA, as necessary. As shown, 
in the interference-limited regime, our algorithm reconfigures 
itself so as to achieve IA whereas, in the low-to-moderate 
interference regime, it chooses interference-myopic MIMO 
transmissions. Furthermore, for all investigated interference 
cases and channel conditions, it was shown that the sum-rate 
of the proposed algorithm is higher than that of all considered 
IA-achieving algorithms. 
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